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Abstract. The half-integer quantum Hall effect in epitaxial graphene is compared 
with high precision to the well known integer effect in a GaAs/AlGaAs heterostructure. 
We find no difference between the quantised resistance values within the relative 
standard uncertainty of our measurement of 8.7 x 10 -11 . The result places new 
tighter limits on any possible correction terms to the simple relation i?x = h/e 2 , 
and also demonstrates that epitaxial graphene samples are suitable for application 
as electrical resistance standards of the highest metrological quality. We discuss the 
characterisation of the graphene sample used in this experiment and present the details 
of the cryogenic current comparator bridge and associated uncertainty budget. 
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1. Introduction 

The discovery 6 years ago of the quantum Hall effect (QHE) in graphene sparked an 
immediate interest in the metrological community [H [2]. The QHE is a fascinating 
macroscopic quantum effect occurring in two-dimensional conductors and relates the 
resistance quantum, h/e 2 only to the fundamental constants of nature, h, the Planck 
constant and, e, the elementary charge [3]. Although the QHE has been used successfully 
in metrology to realise the resistance scale for more than two decades [I], graphene is 
a material with properties like no other. Graphene, a single layer of carbon atoms 
in a hexagonal crystal lattice structure, is a truly two dimensional metal with a 
linear dispersion relationship characteristic of massless Dirac-type charge carriers [5]. 
The unique bandstructure of this semi-metal has both practical and fundamental 
implications. Firstly, the massless nature of the charge carriers leads to a Landau 
level spectrum with an energy gap between the first two levels which is around 5 times 
larger than that in semiconductor materials for magnetic fields around 10 tesla. This 
implies that the QHE in graphene can be observed at much reduced magnetic fields 
and/or much higher temperatures [6]. Secondly, the marked difference in bandstructure 
and charge carrier characteristics between graphene and semiconductor systems allows 
for a demonstration of the universality of the quantum Hall effect through a rigorous 
test of the material independence of the value of Rk, the von Klitzing constant. 

Theory predicts no measurable corrections to the simple relation Rk = h/e 2 . The 
quantum Hall resistance is considered to be a topological invariant [7] , not altered by the 
electron-electron interaction, spin-orbit coupling, or hyperfine interaction with nuclei. 
It has also been shown that the quantised Hall resistance is insensitive to much more 
subtle influences of the gravitational field [8]. Recently, a quantum electrodynamical 
approach to charge carriers in a magnetic field has predicted a tiny correction to the von 
Klitzing constant of the order of 10 -20 for practical magnetic field values [9]. However, 
the size of this predicted correction is about 10 to 8 orders of magnitude smaller than the 
most accurate measurement techniques available and therefore untestable. Nevertheless, 
the fundamental nature of the Hall resistance quantisation makes experimental tests of 
its universality of the utmost importance, in particular, for improving our knowledge 
of two fundamental quantities of nature: the electron charge and the Planck constant. 
The precision obtained through a universality test as presented here is much greater 
than is possible by a comparison of the best values of the constants h, e, and Rk [10J . 
Analysis of the complete set of published results carried out by CODATA [10] showed 
no deviation from h/e 2 to within 2 x 10 -8 , which calls for more accurate measurements. 
Universality of Rk will strongly support the pending redefinition of the Si-units for mass 
and current in terms of h and e |llj . 

A direct comparison of the Hall resistance in two different substances does not 
prove the exactness of the relationship Rk = h/e 2 ; however, material independence 
is a significant factor in establishing the fundamental nature of Rk- This material 
independence turns out to be rather difficult to establish. Indeed the characteristics 
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of QHE samples must satisfy very stringent requirements |T2] and in 30 years 
only silicon MOSFETs (metal-oxide-semiconductor field-effect transistors) and III-V 
(GaAs/AlGaAs or InGaAs/InP) heterostuctures did so. 

The first accurate measurements of the QHE in graphene were performed by 
Giesbers et al. p3] on exfoliated samples. The precision obtained in these measurements 
was 15 ppm and limited by the high kQ) contact resistances together with a 
small (~ 2.5 yuA) maximum source-drain current which these samples could sustain 
before breakdown of the QHE occurred. A large measurement current determines the 
maximum signal-to-noise ratio and increasing this breakdown current is key to high 
resolution measurements. One established method of increasing the breakdown current 
is to increase the sample width j3] which is not easy to achieve with the exfoliation 
technique. 

In 2009 a number of groups around the world almost simultaneously succeeded in 
growing large-area wafers of epitaxial graphene by sublimation of SiC with a quality 
good enough to observe the quantised Hall resistance [HI EEl [H3, HH [18]. In an 
indirect comparison with a GaAs/AlGaAs device via an intermediate room-temperature 
standard resistor, we demonstrated universality of Rk with an uncertainty of 3 parts 
in 10 9 [16]. The measurement system [19] was identical to that used by Giesbers et 
al. [13] and the key factors in the improvement were the very low contact resistances 
and large sample size which were achieved, resulting in an order of magnitude increase 
in the breakdown current. Subsequently we undertook a direct comparison between 
epitaxial graphene and two GaAs/AlGaAs devices using a modified measurement system 
demonstrating equivalence to an accuracy of 8.6 parts in 10 11 [20]. Recently, much 
progress has been made on metrological quantum Hall measurements in exfoliated 
graphene devices with accuracies continuously increasing to around 5 parts in 10 7 
[2H[22]. 

In this paper we present a detailed description of the graphene device 
characterisation, measurement system, data analysis and resulting uncertainty budget 
which underpins the direct comparison result [20] . 

2. Characterisation of epitaxial graphene 

In [16] we reported the first accurate measurements of the quantum Hall effect in large 
high-quality epitaxial samples. The material studied in our experiments was grown on 
the Si-terminated face of a 4H-SiC(0001) substrate [23J. The reaction kinetics on the 
Si-face are slower than on the C-face because of the higher surface energy, which helps 
homogeneous and well-controlled graphene formation. Graphene was grown at 2000 °C 
and 1 atm Ar gas pressure, which result in monolayers of graphene atomically uniform 
over more than 50 /im 2 , as shown by low-energy electron microscopy. Twenty Hall bar 
devices of different sizes, from 160 /xmx 35 /im down to 11.6 /imx 2 /im were produced on 
each 0.5 cm 2 wafer using standard electron beam lithography and oxygen plasma etching 
(Figured]). Atomic force microscopy (AFM) images revealed that the graphene layer 
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Figure 1. SEM image of a typical device and contact labels. The values for the 
parameters L, W, I and l p are listed in Tab. [2j 



covers the substrate steps like a carpet, preserving its structural integrity [16]. Contacts 
to graphene were produced by straightforward deposition of 3 nm of Ti and 100 nm of 
Au through a lithographically defined mask followed by lift-off, with a typical area of 
graphene- metal interface of 10 4 /zm 2 for each contact. Using transport measurements in 
low magnetic fields, we established that the manufactured material was n-doped, with 
the measured electron concentration in the range of (5.5 — 9) x 10 11 cm -2 , mobility about 
2400 cm 2 V _1 s _1 at room temperature and between 4000 and 7500 cm 2 V~ 1 s~ 1 at 4.2 K, 
almost independent of device dimensions and orientation with respect to the substrate 
terraces. All results presented in this work were obtained on the largest Hall bar (i.e. 
160 x 35 H- 

In fact epitaxial graphene grown on SiC is always strongly n-type doped with carrier 
densities typically in the range of 10 12 — 10 13 cm -2 . The doping of the graphene is 
caused by the so-called "dead layer" of carbon atoms in between the SiC substrate and 
graphene. This layer is non-conducting and characterised by a 6\/3 x 6\/3 supercell of the 
reconstructed surface of sublimated SiC. Missing or substituted carbon atoms in various 
positions of such a huge supercell in the dead layer create localised surface states with a 
broad distribution of energies within the bandgap of SiC. It appears that the density of 
such defects is higher in material grown at low temperatures (1200 — 1600°C) resulting 
in graphene doped to a large electron density [21]. On the other hand, the material 
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Figure 2. Transverse (p xy ) (contacts 4 and 5 in Fig. Q} and longitudinal (p X x) 
(contacts 6 and 7) resistivity measurement. In all measurements in this work the 
source-drain current is between contacts and 1. The horizontal lines indicate the 
exact quantum Hall resistivity values for filling factors v = ±2 and ±6. 



used in this work is grown at high temperatures in a highly pressurised atmosphere of 
Ar which seems to improve the integrity of the reconstructed 'dead' layer, leading to a 
lower density of donors on the surface and, therefore, producing graphene with a much 
lower initial doping [TU E5] . 

Figure [2] shows a typical measurement of the longitudinal and transverse resistivity 
as a function of magnetic flux density B at low temperatures. Wide plateaux are 
observed in the transverse resistivity which are accompanied by a vanishingly small value 
of the longitudinal resistivity. The sequence of quantum Hall plateaux and absence of 
the v = 4 plateau confirms that the sample is monolayer graphene. At zero field a small 
weak localisation peak is visible - a signature of quantum coherence in disordered Fermi 
liquid in epitaxial graphene [26]. Charateristic scattering lengths have been determined 
from the analysis of the weak localisation correction to the Drude conductivity. For 
our graphene grown on SiC we find that the phase coherence length L^, w 1 fim at 
low temperatures, the intervalley scattering length Li »s 0.2 fim and the intravalley 
scattering length « 0.02 /xm, the latter two being virtually temperature independent. 
These experimental scattering lengths set a limit of at least 50 ps on the spin relaxation 
time in this material. 

For accurate quantum Hall measurements it is important that the contact resistance 
is well below m 100 Q [I]. Table [j] shows the measured contact resistances for our device 
and demonstrates that this condition is satisfied for most of the contacts. Our graphene 
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Table 1. Three terminal contact resistances measured at B = 14 T in the v = 2 
quantum Hall state. The measurement current was 10 /zA and a lead resistance of 
2.7 SI was subtracted. 
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device was approximately two years old at the time of these measurements and had been 
thermally cycled between 300 mK and room temperature more than a dozen times. Also 
the contacts have been re-bonded several times in the T08-header which has caused 
damage to some of the bonding pads and is the most likely cause of the poor resistance 
of contacts 3 and 5 (on pristine devices we routinely get all contact resistances below 
1 fi). 

In order to extend the range of applicability the v = 2 quantum Hall plateau 
can be brought down in magnetic field by reducing the carrier density using a novel 
photochemical gating technique [25]. This is achieved by coating the device with a 
polymer bilayer, a spacer layer (PMMA) followed by an active polymer (ZEP520) able 
to generate potent electron acceptors under UV light exposure. The effect can be 
reversed by heating the sample to 170°C just above the glass transition temperature 
of the polymer. Using this technique the carrier density of our device was reduced to 
4.6 x 10 11 cm -2 resulting in an extremely wide v = 2 quantum Hall plateau as can be 
seen in figure [3]^b). 

Accurate quantum Hall resistance measurements require that the longitudinal 
voltage remains zero (in practice, below the noise level of the voltmeter) to ensure 
the device is in the non-dissipative state, which can be violated by the breakdown of 
the QHE at high source-drain current levels. Fig. [3]^a) shows the determination of the 
breakdown current I c for different values of B along the v = 2 plateau. Here we define I c 
as the source-drain current, I S( i, at which V xx > 100 nV. Above the breakdown current 
the voltage increases rapidly as a function of current as a consequence of the sudden 
increase in electron temperature and subsequent increase in longitudinal resistivity j27j . 
In figure [3]^b) the values of the breakdown current are plotted as a function of B together 
with p xx and p xy . 

Figure H] shows the measurement of V xx as a function of I S( i at 11 and 14 T in more 
detail (these are the magnetic flux densities at which the comparison measurements 
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Figure 3. (a) Measurement of V xx (contacts 4 and 7) as a function of source-drain 
current at different values of magnetic flux density ranging from B = 7 T to (left 
hand curve) to 14 T (right curve) in steps of 1 T. (b) Transverse (p xy ) (contacts 4 and 
5) and longitudinal (p xx ) (contacts 4 and 7) resistivity measurement at the reduced 
carrier density of4.6x 10 11 cm" 2 measured at I s d = 1 /iA together with the measured 
breakdown current, I c . Dashed blue line indicates position of the exact v = 2 filling 
factor for the low field carrier density. 
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Figure 4. (a) Log- log plot of V xx as a function of source-drain current at 11 T (Green 
squares) and 14 T (Black dots). The red lines are fits to the data at high source-drain 
currents. 5 nV is the noise floor of the voltmeter. 

are performed). On the log- log plot it is clear that above the breakdown current the 
longitudinal voltage shows a power-law dependence on current. The red lines in the 
figure are a fit to this part of the trace. For the highest flux density of 14 T the 
breakdown current is ~ 500 fiA, however, the maximum source-drain current used in 
our comparison measurements is 100 fiA. The extrapolated fitted line indicates that 
the expected V xx is ~ 1.3 pV at 100 fiA and therefore R xx ~ 13 nf2. Note that close 
to the breakdown current V xx drops away much more quickly than this exponential 
behaviour before it disappears in the noise of the DVM (±5 nV) (the shape of the 
breakdown curve is remarkably similar to that observed by Cage et al. [28] for GaAs 
heterostructure devices) and therefore the R xx ~ 13 nf2 is the upper bound. On the 
basis of this analysis we can be confident that the longitudinal resistivity is vanishingly 
small and that there is a considerable safety margin on I S( i in our experiments. A few 
measurements were performed at a lower magnetic flux density around 11 T. For these 
measurements a lower of 60 fiA was chosen and a smilar extrapolation of the high 
source-drain currents gives V xx ~ 20 pV and R xx ~ 0.3 fifl. 

Returning to figure [3]^b), we note that the measured I c for contact pairs on either 
side of the device is virtualy identical. I c as a function of magnetic flux density shows a 
rather unusual dependence as the breakdown current continues to increase in the v = 2 
quantum Hall state reaching « 500 fiA at our maximum field of B = 14 T. This 
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Figure 5. Middle: Measurement of A GaAs/A1GaAs _ Graphonc = 

[i?ff(GaAs/AlGaAs, T = 1.5 K) - i? H (Graphene, T)]/(R K /2) and p xx as a 
function of temperature for the graphene device. B = 14 T for the graphene device 
and 10.5 T for the GaAs/AlGaAs device. The measurement current was 60 /xA. 
Top and bottom: high-resolution measurements of A GaAs / A i GaAs _ Graphene and p xx 
demonstrating ppb-level quantisation. For the middle panel the carrier density was 
4.6 x 10 11 cm" 2 and for the top and bottom panel the density was 6.7 x 10 11 cm~ 2 . 



behaviour is very different from that observed in conventional semiconductor systems 
where the breakdown current peaks at the exact integer filling factor [I] indicated by 
the dashed blue line in Fig. [3]^b). This behaviour can be explained by a magnetic field 
dependent charge transfer mechanism between the interfacial layer and graphene [30] . 
The amount of charge transfer is determined by the interplay of classical (as a result 
of the geometry) and quantum (as a result of the density of states) capacitance of the 
device structure. The density of states oscillates as the Landau levels (LL) pass through 
the Fermi level when the magnetic field increases. The result is that the carrier density 
in the graphene layer is also oscillating and specifically when the Fermi level is between 
the N = and N — 1 LL, the carrier density increases linearly with magnetic field. This 
effectively pins the filling factor at v = 2 over a broad range of magnetic field resulting in 
a novel quantum Hall state which is ideally suited for high precision resistance metrology. 
The anomalous pinning is responsible for the extremely high breakdown current and wide 
operational parameter space of an epitaxial graphene quantum Hall device. 

Figure [5] demonstrates the robustness of the v = 2 quantum Hall state as a function 
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Figure 6. Variation of AcaAs/AiGaAs-Graphene as function of p xx for three different 
charge carrier densities determined from data such as displayed in fig. [5j Red = 
9.1 x 10 11 cm -2 , Green = 6.7 x 10 11 cm -2 and Black = 4.6 x 10 11 cm" 2 . Open black 
triangles are measured in reverse field direction. Solid lines are fits to the data. 

of temperature. Here A Ga As/AiGaAs-Graphcne = [^(GaAs/AlGaAs, T = 1.5 K) - 
Rh (Graphene, T)}/ (R K /2) quantifies the difference between the graphene sample and a 
fixed reference GaAs/AlGaAs device (the measurement technique is explained in detail 
in section H]). In a separate measurement the longitudinal resistivity was measured for 
the graphene device and plotted in the same graph. The temperature for the reference 
device is held constant at 1.5 K while the temperature for the graphene device is varied. 
The middle panel of fig. [5] shows a low resolution measurement across a wide temperature 
range and demonstrates that Rh starts to deviate as soon as a measurable p xx appears. 
The top and bottom panels are high resolution measurements which demonstrate ppb- 
level quantisation up to 15 K when p xx is of the order of several tens of /iV's (the high 
resolution measurements where performed on an earlier cooldown cycle of the device 
when the charge carrier density was higher). 

In figure ED we have plotted A GaA s/AiGaAs-Graphcne versus p xx for three separate 
cooldown cycles of the device. Each time A GaAs/AiGaAs- Graphene is proportional to 
p xx , identical to the well-known empirical relationship ARh = kp xx which has been 
observed for traditional semiconductor systems jl]. The value of k varies significantly 
between different runs which are signified by different charge carrier densities (in between 
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cooldowns the sample was exposed to UV radiation in order to reduce the carrier density 
in small controllable steps). For the final and lowest charge carrier density we also 
measured k in opposite magnetic field direction. The fact that the sign of k changes 
with the 5-field direction indicates that the observed relation between Rh and p xx is 
due to inhomogeneity of the charge density in the sample as proposed by van der Wei 
et al. [29] rather than the finite width of the voltage probes for which no sign change 
is expected. In the latter case k = ^ « 0.7 for both B- field directions which is clearly 
inconsistent with the experimental data. A plausible explanation for the increase in k 
with reducing carrier density could be the increase in inhomogeneity in the sample as 
we get closer to the Dirac point [I] . This is supported by the fact that the value of k was 
found to be different for different combinations of measurement contacts (not shown). 

Combining the results from figure H] where we estimated R xx with the value of 
k obtained in figure [6] allows us to estimate the expected relative error in Rh- For 
B = 14 T and I s d = 100 fiA we obtained R xx = 13 nQ which with k = 0.16 implies 
that we can expect a relative error in Rh much less than 10~ 12 . Similarly for 60 p,A at 
B — 11 T the relative error would be much less than 10~ n . 

In figure [7] (middle) the same experiment as in figure |5]is repeated but this time with 
magnetic flux density as the parameter and T = 300 mK. The magnetic flux density for 
the reference device is held constant at the centre of the v = 2 plateau and the magnetic 
flux density for the graphene sample is varied. Top and bottom panel are high resolution 
measurements which demonstrate ppb-level quantisation (again measured on a different 
cooldown cycle). The graph shows that the v = 2 quantum Hall plateau extends over a 
range of at least 4 T and the measurement is only limited from above by the maximum 
available magnetic field. The inset is a log-log plot of AcaAs/AiGaAs-Graphcne as a function 
of p xx similar to fig. [6] and the blue line is ARh = kp xx with k = 0.39. From the plot it 
can be seen that the data do not quite follow a linear dependence as might be expected. 
A likely reason for this deviation is the fact that as a function of B the charge carrier 
density in epitaxial graphene is not constant which could result in a variation of k as a 
function of B. 

Comparing the performance of our graphene device with that of GaAs devices shows 
a significantly wider operational parameter space for the v = 2 quantum Hall state. In 
GaAs the field range is usually a few tenths of teslas up to 1 Tesla for the best devices 
and for temperature the operational range is generally below 2 K j4J [38]. This result 
is not surprising given the fact that the energy spacing between the n = and n = 1 
Landau levels is approximately seven times larger in graphene than in GaAs around 
a magnetic flux density of 10 T. Also the breakdown current in graphene compares 
favorably with that in GaAs. Although a breakdown current as large as 500 /xA has 
been obtained in mesoscopic GaAs devices, our device is only 35 /im wide, at least 
a factor of ten narrower than standard Hall bars used for metrological measurements. 
This leaves significant scope for further improvement if wider, homogeneous, graphene 
devices can be produced. It is important to note that for the lowest charge carrier 
density (~ 4 x 10" 11 cm -2 ) in our graphene device we have not yet reached the centre 
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Figure 7. Middle: Measurement of -A GaAs/A1GaAs _ Graphene = 

[i? ff (GaAs/AlGaAs,B = 10.5 T) - R H (Graphene, B)]/(R K /2) and p xx as a 
function of B on the graphene device. T = 0.3 K for the graphene device and 1.5 K 
for the GaAs/AlGaAs device. The measurement current was 60 /iA. Top and bottom: 
high-resolution measurements of — AcaAs/AiGaAs-Graphcnc and pxx demonstrating 
ppb-lcvcl quantisation. For the middle panel the carrier density was 4.6 x 10 11 cm -2 
and for the top and bottom panel the density was 6.7 x 10 11 cm -2 . Inset: log-log 
plot of the variation of A GaAs / A i GaA s-Graphene as function of p xx for the data in the 
middle panel. 



of the v = 2 quantum Hall plateau and the operational parameter space is likely to be 
even wider at higher magnetic flux densities or in devices with even lower charge carrier 
density. 

3. Characterisation of GaAs samples 

The GaAs/AlGaAs heterostructures used in this work were standard Hall bar devices 
in use for routine quantum Hall metrology for the last 20 years at NPL and BIPM. One 
device was supplied originally by PTB (marked as PTB2) and the other by Laboratoires 
d'Electronique Philips (LEP). The dimensions and characteristic parameters are listed 
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Table 2. Device parameters. The dimensions are defined in Fig. Q] The contacts on 
the PTB2 sample were made of small tin balls at the edges of the chip and so voltage 
probe width is not applicable. 
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(xlO 11 cm" 2 ) 


PTB2 
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40 


4.6 
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2.2 


0.4 
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50 
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0.16 
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0.06 


0.024 


0.75 


4.6 



in Tab. [Altogether with those of the graphene device discussed in the previous section. 

Before commencing the high-accuracy measurements both GaAs/AlGaAs samples 
were characterised according to the guidelines for quantum Hall resistance metrology [12] 
(i.e. we confirmed that the three-terminal contact resistance measured on the v = 2 
plateau was of the order of a few ohms for all contacts used and that the longitudinal 
resistivity at the measurement current was below 10 The breakdown current, I c , 

measured at 1.5 K was 150 fiA for the PTB2 device and 100 fiA for the LEP device. 

4. Cryogenic current comparator bridge 

4-1. Design principle 

The measurements were made with a cryogenic current comparator (CCC) bridge which 
has been described in detail in Refs. [T9l l3Tj and is illustrated in simplified form in 
Fig. [HI Isolated current sources 1 and 2 separately drive current through samples Si 
and S2 and associated windings A and B on the CCC. The current ratio can be set 
via electronics to a few parts in 10 6 and this ratio is improved to a level of 1 part in 
10 11 by forming a negative feedback loop from the SQUID (Superconducting Quantum 
Interference Device) sensing the net flux in the CCC to one of the current sources. Our 
standard CCC uses the model A20 nanovolt amplifier manufactured by EM Electronics. 
The performance of this amplifier can be characterised as a voltage noise with Allan 
deviation of 0.14 nV and a current noise with Allan deviation of 0.2 pA for a 10 s 
observation (also see Fig. 4 in Ref. [31]). When measuring two quantum Hall devices 
with a combined impedance of ~ 25 kQ this current noise corresponds to 5 nV and 
dominates the resolution of the measurement system. We therefore chose to use a 
second CCC as a null detector [32]. The potential contacts on S\ and S2 are closed in 
a loop via winding C on this second CCC. This device is configured with just a single 
winding to measure a current null rather than two windings to establish a current ratio. 

The CCC for the bridge current ratio is fitted with an RLC filter which is a modified 
version of the RC filter described in [19]. The modification gives a lower bandwidth of 
1 kHz and was found to be necessary due to the larger interference experienced in this 
experiment involving four cryostats and long connecting cables between laboratories. 
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Also, as a much higher accuracy was required for this precision comparison of quantum 
Hall samples than is needed for routine resistance measurement a guard, driven by a 
unity-gain operational amplifier, was added to the capacitors in the filter on each side of 
the bridge [33] to reduce any leakage across the CCC windings to a negligible level. The 
CCC for null detection was also fitted with a filter (at room temperature) to limit the 
bandwidth to approximately 16 kHz. The filter was constructed using 100 Q resistors 
so as to keep the additional Johnson noise in the null detector circuit at the nV level. 
The potential difference between the bridge circuit and the overall measurement system 
screen in the region of the null detector is held close to zero by the 100 kSl resistor 
connected to the screen in the filter for current source 2. Any leakage resistance in the 
filter capacitors will therefore have a negligible effect and a guard circuit for these is not 
necessary. 

The parameters for both CCC's are listed in Tab. EJ The current sources and 
SQUID control units are connected to a computer which supervises the setting of the 
measurement currents and collects data from the two SQUIDs. Data are collected 
alternately in forward and reverse current direction so as to eliminate electrical offsets 
and drift. 

4-2. Bridge sensitivity 

The CCC bridge sensitivity can be derived from the schematic network in Fig. |HJ At 
bridge balance the flux in the ratio CCC is zero so that I\ — — I. A deviation between 
the ratio R1/R2 from unity will result in a small current, Ind m the null detector circuit. 
We have, 



Ri{I-Ind) = R*{I + Ind), (1) 
R2 _ I — Ind ^ ^ 2Ind ^ 
Ri I + Ind I 
The sensitivity of the null detector was checked by placing a 10 GQ shunt resistor 
in parallel with one of the quantum Hall devices. For a measurement current of 100 /iA 
a SQUID signal for the null detector CCC corresponding to Ind = 70 pA was observed 
which, according to equation [2} corresponds to a deviation of the resistance ratio from 
unity of 1.4 x 10~ 6 . This is in reasonable agreement with the calculated deviation of 
1.3 x 10" 6 , confirming the null detector sensitivity. A second feedback loop can be 
operated to bring the SQUID signal from the null detector to zero using a third winding 
on the CCC [19]. However, this is only necessary when an accurate resistance ratio 
measurement for resistors with deviations from nominal is required. In this experiment 
we are comparing two quantum Hall devices which are expected to be exactly equal 
to a very high degree, so this second feedback loop is not necessary and the open loop 
output of the null detector is simply recorded instead. 
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Figure 8. Schematic of the cryogenic current comparator bridge circuit 
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Table 3. Summary of individual noise components in the bridge circuit and total 
expected noise. The noise attributed to the two CCCs is measured at 10 s and converted 
to voltage noise across one QHR device in case of the current CCC and across two QHR 
devices for the voltage CCC. The noise arising from the resistors is calculated. Last 
two columns are the Allan Deviation (AD) at a 10 second measurement time. Note 
that for white noise of power spectral density ho V 2 Hz _1 the corresponding Allan 



deviation 


is (^ /2r)- 1 /2. 








Component 


Comment 


Sensitivity 


CCC null 


A20 null 








AD (nV) 


AD (nV) 


Ratio CCC 


1600:1600 t 


16 /xA-turn0Q 1 


1.3 


1.3 


Null CCC 


2500 t 


13 /xA-turn^ 1 


1.3 




A20 null 


Voltage noise 






0.14 


A20 null 


Current noise 






5.0 


QHR1 


300 mK 




0.1 


0.1 


QHR2 


1.2 K 




0.2 


0.2 


Filter 


200 n 




0.4 




Total 






1.9 


5.2 



4-3. Noise and measurement resolution 

The ultimate resolution of the bridge can be estimated by considering the noise 
components due to the current balance CCC, the null detector CCC and filter network, 
and the Johnson noise from the quantum Hall devices being measured. Here we use the 
Allan deviation of the time domain signals [34"| [T9] to analyse the various components. 
The design of the bridge is based on a noise performance of the DC SQUID and CCC 
combination with an Allan deviation of 10 /x0o for a 10 s measurement interval. Table [3] 
lists the parameters for the different components of the bridge and the calculated voltage 
noise in the null detector loop in terms of the Allan deviation for a 10 s measurement 
interval. The total noise voltage for the bridge with CCC nanovoltmeter is 1.9 nV 
compared to 5.2 nV when using the A20 nanovoltmeter. 

Figure [9] shows the measured Allan deviation of the CCC null detector with two 
QHR devices connected to the bridge (the master and slave current sources and current 
balance CCC are disconnected). It is again expressed as a voltage noise in the null 
detector circuit. It can be seen from this graph that at 10 s the nanovoltmeter noise is 
approximately 8 nV compared to design value of 1.3 nV. For comparison the measured 
Allan deviation is also shown for the A20 null detector which is similarly ~8 nV at 10 s 
rather than the expected 5 nV. This suggests that both the CCC null detector and the 
A20 instrument are experiencing a similar level of excess noise. The enhancement of the 
noise is not surprising given the fact that we are using two quantum Hall systems each 
with their own magnet power supply and peripheral pumping equipment. In addition the 
connecting cables between these systems and the measurement bridge span a distance 
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Figure 9. Allan deviation of the CCC nanovolt meter (orange) and A20 nanovoltmeter 
(blue) when two QHR devices are connected to the bridge. Also shown is SQUID noise 
of the current balance CCC which was scaled to 1600 turns. 



of more than 10 meters across two laboratories. 

In Fig. [9] the noise from the current balance CCC is also shown and has been scaled 
for a 1600 turn winding and a single quantum Hall resistance to give an equivalent 
voltage noise in the null detector circuit. The Allan deviation goes through a minimum 
of ~ 3nV at about 10 s which is again more than the expected 1.3 nV and then increases 
at one decade per decade in sampling interval. (Note that the minimum in the Allan 
deviation determines the optimum data acquisition time of ~ 10 s in our measurement 
system.) Both SQUIDs exhibit a similar behaviour for time intervals longer than 
30 s. The excess noise in the null detector circuit is clearly the limiting factor in our 
measurement system at present. Excess noise has also been seen in measurement systems 
of this type [32] and whereas in a manual data acquisition system the user can select 
apparently lower noise results, here our system runs for many hours unattended with 
the consequence that all data points are included in our analysis. From Fig. [9] we can 
conclude that the optimal measurement time is approximately 20 s for a given current 
direction. For longer times both the CCC null detector noise and current balance CCC 
noise contributions will start to increase with the null detector at present dominating 
until time intervals longer than 50 s. A future improvement in this experiment would 
be to select a graphene and GaAs device with coinciding v = 2 plateaux which can be 
co-located in the same cryostat and thereby reduce the complexity of the connecting 
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Figure 10. Typical Allan deviation of the resistance ratio for a 1:1 measurement of 
graphene against GaAs/AlGaAs heterostructure at v = 2 for a measurement current 
of 100 fxA. 



cables. 

Figure [10] shows the Allan deviation for a 3.5 h measurement of graphene against 
GaAs at v — 2. The measurement data have been analysed in blocks of three intervals 
consisting of a 1 s settling time plus 10 s data collection interval (for optimum CCC 
performance) and a 1 s data transfer each to give a value for the resistance ratio every 
36 s. (Note that the settling time for the CCC nanovoltmeter is much faster than that 
of the A20 where we have to allow for at least 10 s settling time [19]. We regularly 
checked the stability of the servo loop by varying the wait-time after current reversal 
and inspecting the measurement result for systematic correlations.) Adding all the 
noise components in Tab. [3] in quadrature and using the observed noise for the CCC 
nanovoltmeter (which dominates) and current balance CCC in Fig. Ogives a total noise 
of 8 nV for a 10 s measurement interval. For a measurement current of ±100 fiA in 
Rk,u=2 this corresponds to a relative deviation of 3 parts in 10 9 . This number can be 
related to a block of three intervals by multiplying by a scaling factor of J3/2 [19] to 
give 4 parts in 10 9 for the expected Allan deviation. From Fig. [10] we see that the first 
data point is also at 4 parts in 10 9 confirming that the overall measurement resolution is 
consistent with the measured noise. The Allan deviation decreases as l/y/r as expected 
for white noise and shows that the current reversal technique efficiently removes the 
non-white noise visible in Fig. [9] After 3.5 h measurement time a relative uncertainty 
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of 2 parts in 10 10 is achieved. If the two CCC's can be made to operate at their 
optimum noise performance, a single block of three measurements would give a relative 
uncertainty slightly better than 1 part in 10 9 and 6 parts in 10 11 after 3.5 h (blue dot 
and dashed line in Fig. [TO]) . 

5. Uncertainty analysis 

5.1. Type B analysis 

5.1.1. Ratio error The ratio error is determined by connecting both 1600 turn windings 
in series opposition. A large (10 mA) current is passed through the windings and the 
SQUID signal recorded (total applied current linkage is 32 A-turn). This process is 
repeated several hundred times in both forward and reverse current direction to improve 
the measurement resolution. A least squares fit to the data gives 17 ± 7 /iV for the 
residual SQUID signal. The sensitivities of the SQUID and CCC are 0.7 V0q 1 and 
16 yuA-turn</>o 1 , respectively, giving an imbalance of 3.6±0.9x 10~ 10 A-turn. In principle 
it is possible to correct the measurements results for this imbalance because it is a fixed 
property of the CCC. However, given the small size we have not done this and have used 
the imbalance as an uncertainty component, resulting in a relative standard uncertainty 
of 1.1 x 10~ n in the uncertainty budget. 

5.1.2. Leakage Leakage resistances must be controlled to a very high level in precision 
resistance measurements depending on where they appear in the bridge circuit. The 
most stringent condition occurs for the wiring connecting the device in the cryogenic 
system to the bridge. A leakage between the high and low potential side of the device 
would appear as a resistance in parallel with the quantum Hall device and result in a 
direct error. Given that the quantum Hall resistance at v = 2 is approximately 10 4 f2, 
this implies a leakage > 10 15 Q for a relative error less than 10 -11 . Measuring such a 
large resistance on long leads (there is no screen between high and low potential leads 
for ~ 3 m from the top of the cryostat probe to the sample in the 14 T/300 mK system 
at NPL) is not straightforward. Using a Keithley 6430 electrometer and a repeated 
±100 V excitation during a 10 hour measurement resulted in 0.1 ± 0.06 pA, giving a 
relative standard uncertainty of 1 x 10~ n in the uncertainty budget. 

5.1.3. Servo error Any residual error in the current ratio between master and slave 
will translate directly into a false reading on the null detector. In order to test and 
correct for this error, we record the residual SQUID voltage from the current balance 
CCC during measurements. From this we can calculate the resulting current in the null 
detector CCC via, 
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and apply a correction. Typically we find corrections in the range of ~ 3 ± 0.6 fA 
which corresponds to a relative correction of 3 ± 0.6 parts in 10 11 to the resistance ratio 
for a measurement current of 100 fiA. The reason this small error arises is that the 
slave current source drifts very slightly relative to the master source after each current 
reversal, due to self-heating in the current source components. The feedback servo has 
a single integral term so a linear drift in the current source translates into a steady error 
voltage at the output of the SQUID. For our system, the residual error current in the 
CCC of ~ 3 fA quoted above corresponds to an open loop drift in the current sources 
of approximately 1 part in 10 6 over the period of a 20 s measurement. 

5.1.4- Miscellaneous In the comparison measurements are performed in sets of ten 
intervals of successive forward and reverse current (rather than the shorter blocks of 
three intervals which were used in the previous section to analyse the time evolution of 
the Allen deviation). Each data set is written to disk and a new measurement started. 
Typically this operates for many hours and hundreds of data sets are collected. Each 
data set of ten intervals is analysed using a least squares fitting routine to give a value 
for the resistance ratio and standard deviation. The residuals of the fit are analysed in 
order to spot potential problems with the system such as excessive settling time, noise 
and interference. 

Figure [TlT a) shows the average of 300 sets of fit residuals (i.e. 10 hours total 
measurement time) and a marked oscillation becomes visible in the data (indicated by 
the green curve). The oscillation starts at the first forward current measurement interval 
and damps out with a time constant of about 60 s. We speculate that this effect is caused 
by a slight heating of the CCC helium dewar due to the measurement current. Normally 
we start a measurement set from zero current and ramp to maximum forward current, 
collect data, ramp to maximum reverse current, etc [blue curve in Fig. fTTT a)]. A current 
reversal takes approximately 1 s. At the end of a measurement set the current is ramped 
back to zero, data written to file and a new measurement started which combined takes 
about 6 s, i.e. the current is zero for 6 s between measurement sets. Because the down- 
leads to the CCC are not superconducting, dissipation occurs during the measurement, 
and drops for 6 s between measurement sets, causing a temperature-pressure oscillation 
in the dewar. Averaging the residuals in Fig. [TTT a) results in an error of ~ 1 x 10~ 10 
which is not insignificant for our comparison. 

Figure [TTT b) displays the average of a similar 10 hour measurement where the 
measurement current was always present in the CCC (i.e. the measurement current was 
left at maximum reverse current at the end of a measurement set and the new set would 
start by sweeping from reverse to forward maximum current). As can be seen from this 
graph the oscillation has disappeared and the mean of all residuals averages to 3 x 10~ 15 
and any remaining error is therefore negligible. 

This effect demonstrates the care that must be taken when making very high 
precision measurements and the power of automated measurements and data analysis 
in tracking down extremely small systematic errors. 
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Figure 11. (a) Average of fitting residuals for 300 measurement sets using our 
standard measurement sequence (blue curve) and (b) modified sequence. Green curve 
in (a) is guide to the eye. 



6. Type A analysis 

The automated measurement system allows us to make continuous measurements over 
many hours (usually during the night) so that the statistics on the final answer can be 
improved. Figure [12] shows a histogram of all CCC null detector values, expressed 
in terms of the voltage at the SQUID output, obtained from individual data sets 
of 10 current reversals. The SQUID output voltage is the correct parameter to 
combine to give a distribution, since the noise recorded by this detector does not 
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Figure 12. Histogram of null detector signal and Gaussian fit. 



depend on the measurement current in the bridge (the data consists of sets at several 
different measurement current values). The total number of data sets for the complete 
measurement campaign was 2300, equating to almost 100 hours of data collection. The 
solid line is a least-squares Gaussian fit to the histogram and gives a value of 245 /iV 
for the standard deviation. Using the SQUID amplifier gain and the sensitivity values 
in Table El the standard deviation of 245 fiV translates to a current Ind of 150 fA. 
Finally, using equation [2] and a measurement current of 100 /iA, this standard deviation 
is equivalent to a relative deviation of 3.1 x 10~ 9 in the resistance ratio. Since the 
measurements were done at a range of currents, an effective number of measurements 
has to be calculated from 



N eff = Znj'/I? 



(4) 



where Jo = 100 /iA, and ni is the number of measurements taken at measurement 
current l{. This equation gives n e // = 1057 in our case and dividing the relative 
uncertainty by the square root of this number leads to a final type A uncertainty for the 
entire data set of 9.4 x 10 -11 compared with 8.6 x 10 -11 obtained in the next section. 



7. Comparison result 



In order to test the universality of the quantum Hall effect, two different devices need 
to be set up at the same quantum Hall plateau so that a one-to-one comparison of 
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Figure 13. Measurement of AcaAs/AiGaAs-Graphene for the direct comparison of 
Rh (GaAs/AlGaAs) and Rh (Graphene) as a function of Isd- The uncertainty bars 
represent the ±lcr standard deviation of the mean. Different symbols are explained in 
the text. The red line is the weighted mean of all the data points and the green lines 
signify ±ler. 



Precision comparison of the quantum Hall effect in graphene and gallium arsenide 24 

resistance can be made (in principle one could also compare different index plateaux if 
different ratios are available on the CCC). Unfortunately, the quantum Hall plateaux of 
the graphene and the GaAs/AlGaAs devices do not overlap and two cryostats have to 
be used for this experiment, one of which was the traveling quantum Hall system of the 
BIPM. Hartland et al. [35] used a clever method of tilting one device with respect to 
the magnetic field direction to make the centers of the plateaux in GaAs and Si occur at 
the same magnetic field, a trick not currently possible in our quantum Hall probe. As 
explained in the previous section using two cryostats does lead to a significant increase 
in noise and subsequent lengthening of the measurement time. 

The graphene sample was mounted in a 14 T/300 mK cryostat and connected to 
the Slave side of the CCC bridge. Two GaAs/AlGaAs samples were mounted in the 
transportable 11.7 T/1.5 K cryostat provided by the BIPM and connected to the Master 
side. Prior to making comparison measurements all samples were fully characterised |12j . 
For the graphene sample the B v=2 was set to 14 T, our maximum available magnetic 
field, and either 9.5 T or 10.5 T for samples PTB2 and LEP, respectively. The winding 
ratio on the bridge was always 1600:1600 turns. The red triangles in Fig. [13] are the 
results for GR9 against PTB2 for 4 different source-drain currents in the devices (using 
contacts 2 and 7 on the graphene device for Rh and contacts 4 and 5 on the GaAs 
devicdjl). The pink diamond is a measurement for GR9 against LEP at at I s d = 50 /iA 
(using the same contact configuration as for the red triangle measurements). Here 
A G aAs/AiGaAs-Graphcne = [Rh (GaAs/ AlGaAs) - R H (Graphene)] / (R K /2) and each data 
point consists of an average of between 3 and 10 hours worth of data. The uncertainty 
increases for lower I s d because the signal-to-noise ratio is worse for lower I s ^. 

To eliminate the possibility of errors due to non-zero p xx we repeated a number 
of these measurements for non-opposite contacts. Note that it is very difficult to 
measure p xx directly to the required level of precision (see for example fig. H]). The 
green dot is a measurement using contacts 3 and 4 for i?#(GaAs/AlGaAs) on the PTB2 
device at I S( i = 100 /iA. The blue square is a measurement using contacts 2 and 4 
for Rh (Graphene) on the graphene device at I s d = 75 /iA. Another test to check for 
small errors is to reverse the direction of magnetic field on the graphene sample. The 
results of this measurement is represented by the light blue hexagon for I s d = 60 /iA 
(using the same contact configuration as for the red triangle measurements). Finally, 
the devices were exchanged between the NPL and BIPM cryostats in order to check for 
small parasitic leakages (the magnetic field was set to 11.5 T for the graphene device in 
the BIPM cryostat). This data point is shown as the black square in Fig. [13] (using the 
same contact configuration as for the red triangle measurements). 

The data in Fig. [13] show no sign of any systematic errors in the measurement 
campaign and so all results can be combined to give a weighted mean of 
A G aAs/AiGaAs-Graphene = (-4.7 ± 8.6) x 10~ n . Combining this with the total type B 
uncertainty determined in Sec. 15.11 gives 8.7 x 10~ n as a final combined uncertainty on 

| The contact labels for the GaAs devices are the same as for the graphene device depicted in fig. [1] 



Precision comparison of the quantum Hall effect in graphene and gallium arsenide 25 
the difference between i?^(GaAs/A!GaAs) and Rh (Graphene 



8. Conclusion 

Previously our knowledge of the universality of the QHE has been limited to the level 
of 2 or 3 x 10~ 10 for comparisons between GaAs and Si or between identical GaAs 
devices [351 ESI [371 EH]- However both GaAs and Si are traditional semiconductors 
with a parabolic bandstructure and governed by the same physics. Graphene is a 
semi metal with a linear bandstructure and is described by Dirac-type massless charge 
carriers and so universality in terms of material independence goes well beyond the 
comparison between two semiconductors. It does directly support the Thouless-Laughlin 
argument [7J that the Hall conductivity is a topological invariant and is a fundamental 
test of condensed matter theory. 

Our results on material independence is the strongest evidence yet that the 
hypothesis that the resistance is quantised in units of h/e 2 is correct and thereby 
supports the pending redefinition of the SI units kilogram and ampere in terms of h 
and e. Note that the correctness of this equation can only be shown in a comparison 
of -Rh with an independent realisation of h/e 2 , for example via the Thompson-Lampard 
capacitance [39]. However, the relative uncertainty which can be obtained in such 
an experiment is some 3 orders of magnitude worse than in a test of the material 
independence. 

In our universality experiment the maximum source-drain current that the GaAs 
device can sustain without dissipation limits the measurement uncertainty, whereas a 
potentially lower uncertainty can be obtained in a consistency check of two graphene 
devices. Also the measurement system could still be significantly improved by reducing 
the excess noise or by adopting a different measurement technique such as demonstrated 
in Ref. [37]. 

The journey from the original discovery of the QHE in graphene to a quantum 
resistance standard which outperforms the established technology in many aspects has 
been remarkably short. For epitaxial graphene the robustness of the quantisation in 
terms of temperature, magnetic field and source-drain current is exceptional. The 
material is cheap and relatively easy to fabricate and process. It allows for the realisation 
of a quantum resistance standard with modest means, e.g. a small superconducting 
magnet and cryocooler. As such it will improve the proliferation of quantum standards 
and allow many smaller laboratories to realise their own resistance scale. One even could 
envisage university students being able to perform QHE experiments, much in the same 
way as the discovery of high-T c superconductors enabled table-top experiments with 
Josephson junctions in many science classes. 

§ This relative uncertainty is fractionally larger (0.1 x 10 -11 ) than the one published in Ref. [20] due 
to a more thorough determination of the systematic components. 
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